Screening in gated bilayer graphene 
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The tight-binding model of a graphene bilayer is used to find the gap between the conduction and 
valence bands, as a function of both the gate voltage and as the doping by donors or acceptors. The 
total Hartree energy is minimized and the equation for the gap is obtained. This equation for the 
ratio of the gap to the chemical potential is determined only by the screening constant. Thus the 
gap is strictly proportional to the gate voltage or the carrier concentration in the absence of donors 
or acceptors. In the opposite case, where the donors or acceptors are present, the gap demonstrates 
the asymmetrical behavior on the electron and hole sides of the gate bias. 

PACS numbers: 73.20.At, 73.21.Ac, 73.43.-f, 81.05.Uw 



Bilayer graphene has attracted much interest partly 
due to the opening of a tunable gap in its electronic 
spectrum by an external electrostatic field. Such a phe- 
nomenon was predicted by McCann and Fal'ko [l( and 
can be observed in optical studies controlled by apply- 
ing a gate bias 0, & II S, @, • In Refs - B Hl> within 
the self-consistent Hartree approximation, the gap was 
derived as a near-linear function of the carrier concen- 
tration injected in the bilayer by the gate bias. Recently, 
this problem was numerically considered flfjl ] using the 
density functional theory (DFT) and including the ex- 
ternal charge doping involved with impurities. The DFT 
calculation gives the gap which is roughly half of the gap 
obtained in the Hartree approximation. This disagree- 
ment was explained in Ref. fdcj a s a result of both the 
inter- and intralayer correlations. 

In this Brief Report, we study this problem within the 
same Hartree approximation as in Refs. [H, 0], but in- 
cluding the external doping. We consider the case, where 
the carrier concentration in the bilayer is less than 10 13 
cm~ 2 , calculating the carrier concentration on both lay- 
ers. Then, we minimize the total energy of the system 
and find self-consistently both the chemical potential and 
the gap induced by the gate bias. Our results completely 
differ from those obtained in Refs. [H, Q even for the 
range where the external doping is negligible. The de- 
pendence of the gap on the carrier concentration, i. e., on 
the gate voltage, exhibits an asymmetry at the electron 
and hole sides of the gate bias. 

The graphene bilayer lattice is shown in Fig. [TJ Atoms 
in one layer, i. e., a and b in the unit cell, are connected 
by solid lines, and in the other layer, e.g., a\ and b±, 
by the dashed lines. The atom a (ai) differs from b (bi) 
because it has a neighbor just below in the adjacent layer, 
whereas the atom b (bi) does not. 

Let us recall the main results of the Slonchewski- 
Weiss-McClure model UJ, |l2j. In the tight-binding 



the form 

^ = ^E eika ^o( a J- r ) 

3 

^ = ^E eikaj ^( a J+ a - r ) C 1 ) 

j 

^ al = -/^J2 e ika ^ (aj + c - r) 

j 

^ bl = ^Yl eJka ^o( a J + c + a - r), 

j 

where the sums are taken over the lattice vectors a.j and 
N is the number of unit cells. Vectors a and c connect 
the nearest atoms in the layer and in the neighbor layers, 
correspondingly. 




model, the Bloch functions of the bilayer are written in 



FIG. 1: Bilayer lattice 

For the nearest neighbors, the effective Hamiltonian in 
the space of the functions ffl contains the hopping inte- 
grals 70, 71,73,74, and A [13|. The largest of them, 70, 
determines the band dispersion near the K point in the 
Brillouin zone with a velocity parameter v. The param- 
eters 73 and 74 giving a correction to the dispersion are 
less than 70 by a factor of 10 (see Refs. [4j, |5|]). The 
parameters 71 and A result in the displacements of the 
levels at K, but A is much less than 71. Besides, there 
is the parameter U induced by the gate and describing 
the asymmetry of two layers in the external electrostatic 
field. This parameter simply presents the potential en- 
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ergy between two layers, 2U = —edE, where d is the in- 
terlayer distance and E is the electric field induced both 
by the gate voltage and the external charge dopants in 
the bilayer. In the simplest case, the effective Hamilto- 
nian can be written as 
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where the matrix elements are expanded in the momen- 
tum k± — ^fik x — k y near the K points. 
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FIG. 3: Electrostatic model; d is the interlayer distance, d u 
is the wafer thickness. 




FIG. 2: Band structure of bilayer 
The Hamiltonian gives four energy bands: 

e 1A ( q )^±U± + U 2 + q 2 + W\ . 

e 2 . 3 ( q )=±tlL + U 2 + q 2 -W 
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and we denote q 2 = (vk) 2 . 

The band structure is shown in Fig. 
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imal value of the upper energy E\ is \J\J 2 + 7 2 
£ 2 band takes the maximal value \U\ at k = and the 
minimal value U = 7i|J7|/V7i + 4C/2 at 9 2 = 2 ^ 2 (7i + 
2U 2 )/(y 2 + 4J7 2 ). Because the value of U is much less 
than 7i, the distinction between U and U is small and 
the gap between the bands £2 and £3 takes the value 2\U\. 

The eigenfunctions C of the Hamiltonian ^ have the 
form 
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where the C norm squared is 
C 2 = [(e n + U) 2 ^ 2 
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The probability to find an electron, for instance, on the 
first layer is \C\\ 2 + | C2 1 2 , as seen from Eqs. (HJ. 

We assume, that carriers occupy only the bands £2,3, 
so the chemical potential /j, and the gap 2|J7| are less 
than the distance between the bands £1 and £ 2 , i.e., 
(|/x|,2|£7|) < 71. The electron dispersion for the £2,3 
bands can be expanded in powers of q 2 : 
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Then, for q 2 ^ 4[/ 2 , we can use the simple relations: 
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Within such the approximation, many observable effects 
can be analytically evaluated for the intermediate car- 
rier concentration, AU 2 71 \J /j, 2 — U 2 7 2 , where we 
neglect the effect of the "mexican hat". 

At zero temperature, for the total carrier concentration 
n and the carrier concentrations n\,% on the layers, we 
obtain 
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where n — y 2 /Trh 2 v 2 = 1.2 x 10 13 cm -2 and x = [J./U. 

In order to find the chemical potential /j, and the gap 
2\U\ at the given gate voltage, we minimize the total 
energy containing both the energy V of the carriers and 
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the energy V/ of the electrostatic field. Instead of the 
chemical potential, it is convenient to use the variable x 
along with U. Electrons in the e 2 band or holes in the 
£3 band contribute in the total energy of the system the 
energy 
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z\Jx 2 - 1 + In (x + \Jx 2 - 1)] . 



The energy of the electrostatic field, 
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V f = ^-{dE 2 + e w d w E 2 J 
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can be written in terms of the carrier concentrations with 
the help of relations (see Fig. [3]): 

4ne(ni — 714/2) — E and4-7re(n — rid) = £wE w , (11) 

where e w is the dielectric constant of the wafer, the neg- 
ative (positive) rid is the acceptor (donor) concentration 
and we suppose that the donors/acceptors are equally 
divided between two layers. 

We seek the minimum of the total energy as a function 
of two variables, U and x, under the gate bias constraint 
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Excluding the Lagrange multiplier and assuming the in- 
terlayer distance to be much less than the thickness d w 
of the dielectric wafer, we obtain the following equation: 
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Let us emphasize, that this equation is invariant only 
with respect to the simultaneous sign change in rii^ and 
rid, that expresses the charge invariance of the problem. 
At the fixed sign of the external doping rid, the gap on the 
electron and hole sides of the gate bias is not symmetrical. 

The derivatives in Eq. (fl"3"|) are calculated with the 
help of Eqs. (J7J>-([TTJ). As a result, Eq. becomes 
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with the function f(x) = In (x + \J x 2 — 1) and the di- 
mensionless screening constant 
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For the parameters of graphene d = 3.35 A 
and v = 10 s cm/s, we get A = 0.447. 

First, consider an ideal undoped bilayer with rid = 0, 
namely, \jirid/Uno\ <ti 1. We obtain a solution, as xq = 
6.2784, only for one sign in Eq. ([14]) determining the 



polarity of the layers [see Eq. JSJ)]. This value gives 
2\U//j,\ = 2/xo = 0.3186 for the ratio of the gap to the 
chemical potential. According to Eq. ((7]), the gap as a 
function of the carrier concentration takes a very simple 
form: 



2\U/n\ = 



271 



1.08 x 10" n meV-cm 2 , (15) 



where the right-hand side does not depend on the gate 
bias at all, but only on the screening constant A. We can 
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FIG. 4: The gap in units of 71 = 0.4 eV versus the carrier 
concentration for the hole doping with concentration rid = 
— 5x 10 12 cm -2 ; the positive (negative) values of n correspond 
to the electron (hole) conductivity. 
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compare Eq. (|15p with the corresponding result of Ref. 
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Both equations are in the numerical agreement at \n\ — 
O.lno 



10 1 



cm 



However, contrary to Eq. JT5J), 



Eq. (|16p contains the carrier concentration in the right- 
hand side giving rise to the more rapid increase in the 
gap with \n\ uq. This increase also contradicts to the 
DFT calculations [lj|. 

For the bilayer with acceptors or donors, rid 7^ 0, Eq. 
(1141) presents a solution w = Un^/^irid as a function of 
x. We obtain, evidently, the large values of w for x close 
to xq = 6.28. In this region of the relatively large \U\, we 
find again with the help of Eqs. ([7]) and (fT4|) the linear 
dependence 
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For the small gap, \Uno/"/irid\ < 1, we obtain differ- 
ent results for the electron and hole types of conductiv- 
ity. For instance, if the bilayer contains acceptors (see 
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Fig. [4} with concentration rid, the gap decreases linearly 
with the hole concentration and vanishes, when the gate 
bias is not applied and the hole concentration equals nd- 
Starting from this point, the gap increases and, there- 
after, becomes again small (equal to zero in Fig. at 
the carrier concentration corresponding to the minimal 
value of the dc conductivity. Therefore, the difference 
observed in Ref. [IH between these two values of car- 
rier concentrations, at the zero bias and at the minimal 
conductivity, gives directly the donor/acceptor concen- 
tration in the bilayer. Then, for the gate bias applied in 
order to increase the electron concentration, the gap is 
rapidly opening with the electron appearance. 

We see, that the asymmetry arises between the electron 
and hole sides of the gate bias. This asymmetry can 
simulate a result of the hopping integral A in the electron 
spectrum 14J. In order to obtain the gap dependence 
for the case of electron doping, rid > 0, the reflection 
transformation n — > — n has to be made in Fig. 0J 

The gap in the vicinity of the minimal conductivity 
point reaches indeed a finite value due to several reasons. 
One of them is the form of the "mexican hat" shown 
in Fig. [2] Second, the trigonal warping is substantial at 
low carrier concentrations. Finally, the graphene electron 
spectrum is unstable with respect to the Coulomb inter- 
action at the low momentum values. For the graphene 
monolayer as shown in Ref. [l6j], the logarithmic cor- 
rections appear at the small momentum. In the case of 
the bilayer, the electron self-energy contains the linear 
corrections, as can be found using the perturbation the- 
ory. The similar linear terms resulting in a nematic order 
were also obtained in the framework of the renormaliza- 
tion group [17 1. 

In conclusions, the gap 2U opening in the gated 
graphene bilayer has an intriguing behavior as a function 
of carrier concentration. In the presence of the external 
doping charge, i. e. donors or acceptors, this function is 
asymmetric on the hole and electron sides of the gate bias 
and it is linear only for the large gate bias. 
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